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We investigate the ground state of one-dimensional few-atom Bose-Bose mixtures under harmonic confine- 
ment throughout the crossover from weak to strong inter-species attraction. The calculations are based on the 
numerically exact multi-configurational time-dependent Hartree method. For repulsive components we detail 
the condition for the formation of a molecular Tonks-Girardeau gas in the regime of intermediate inter-species 
interactions, and the formation of a molecular condensate for stronger coupling. Beyond a critical inter-species 
attraction, the system collapses to an overall bound state. Different pathways emerge for unequal particle num- 
bers and intra-species interactions. In particular, for mixtures with one attractive component, this species can be 
viewed as an effective potential dimple in the trap center for the other, repulsive component. 

PACS numbers: 03.75.Mn, 67.60.Bc, 67.85.-d 



I. INTRODUCTION 

Cold atoms have become an important tool to create and 
study strongly correlated quantum systems QUI]. One main 
reason is that it is possible to experimentally tune the effec- 
tive low-energy interaction strength between the atoms using 
Feshbach resonances H. This has proven useful particularly 
for Fermi gases 0], whereas for bosons the creation of strong 
interactions is limited by three-body collisions. However, in 
lower (here: one) dimensions there are also other possibilities 
of achieving effectively strong correlations - e.g., by lower- 
ing the atom-number density J5[] and via confinement-induced 
resonances, which exploit the parametric dependence on the 
transverse trapping potential [6]. This allows one to prac- 
tically adjust the coupling strength all the way from infinite 
attraction to hard-core repulsion. 

For single-component bosons in one dimension (ID), the 
extreme case of infinitely repulsive interactions is known as 
the Tonks-Girardeau gas, which has been realized experi- 
mentally 121 |8fl. Here the system maps to an ideal gas of 
fermions, in the sense that the exclusion principle emulates 
the effect of hard-core repulsion ]9|]. The microscopic mech- 
anism of the crossover from the weakly interacting Bose gas 
to the above fermionization limit has been investigated in de- 
tail E3 [Hill H. By contrast, theground state 
for strong attraction is an iV-atom molecule (lql . However, 
exotic fermionized excitations exist for sufficiently attractive 
interactions lfl7i ITsl [Tl . 

In the case of two (or more) bosonic components, a plethora 
of configurations exists: On top of varying both intra- and 
inter-species interactions, also the trapping potentials may 
be made species dependent. Moreover, the experimental 
availability of different two-component mixtures (involving 
not only different hyperfine components Il20l 12111 or isotopes 
122, 12311 . but altogether different atomic species like K-Rb 
024ll ) adds another degree of flexibility. For two ID Bose 
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gases with inter-species repulsion, a generalized, compos- 
ite fermionization exists which may lead to demixing of the 
two components atom by atom J25l l26l l27l l28tl . In a lattice 
pote ntial, even more complex patterns have been found, cf. 
lH 53 ill 51 andRefs. therein. 

In this work, we are interested in the binding between two 
bosonic species, i.e., the crossover from weak to strong inter- 
species attraction. Here little is known except for a gen- 
eral classification based on the harmonic-fluid approximation 
iSH. For fermions, pairing between the two components has 
been predicted, which then form a Tonks-Girardeau gas of 
molecules [34]. Similarly, pairing has been found in attractive 
Bose-Fermi mixtures in various settings SU, 53 IH . While, 
in principle, a ID fermionic component maps to a strictly 
fermionized bosonic one, the physics of realistic Bose-Bose 
mixtures differs in two ways: For one thing, the finite intra- 
species repulsion must compete with strong inter-species at- 
traction. More generally, in contrast to fermions all possible 
/nfra-species interactions are possible and make for interest- 
ing phases. The key goal of this work is to demonstrate effects 
due to the interplay between intra-species and inter-species 
forces. 

Our paper is organized as follows. Section HH introduces 
the model and briefly reviews the computational method. 
The pairing between repulsive components is elucidated in 
Sec. Hill first for the case of a mixture of balanced compo- 
nents (Sees. IIII Al IIIIBb , complemented by a discussion of 
atom-number imbalances and unequal intra-species interac- 
tions (Sec. IIII Cb . Section |IV] deals with the question of how 
the presence of attractive components alters the picture. 



II. MODEL AND COMPUTATIONAL METHOD 

Model The object of investigation is a two-component 
Bose gas (denoted a G {A, B}) subjected to a one- 
dimensional confinement, where the external potentials act- 
ing on the different species are assumed to be the same, i.e. 
U a {x) — U (x). The two species can be considered as two 
internal states (pseudospin |f) and of the same kind of 
Bose atoms, or as different isotopes with the mass m a sa m. 
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In the subsequent sections we denote the atom number of each 
species with N a and the total number with N = Na + Ng. 
This kinematically one-dimensional system of trapped bosons 
can be described in the low-energy limit by an effective one- 
dimensional Hamiltonian with contact interactions. The sec- 
ond quantized Hamiltonian H then reads 



h = fdx E {*£( 



1 d 1 

2 Tjx 1 



U(x) ^ a {x) 



+ f *t {x) & a ( x )* a {x )* a 

+9AB J dx& A (x) ¥ B (x) * B (a;) 9 A (x) , 

where the field operator ^> a (x) (^ (x)) annihilates (creates) 
a boson of the a-species at the position x. The effective intra- 
and inter-species couplings g a and g A B characterize the in- 
teraction between the atoms and can be controlled experimen- 



tally by the scattering lengths and a^ 1 "' , respectively, in 
analogy to the single species case (01 • Furthermore, the stan- 
dard rescaling procedure to harmonic oscillator units has been 
carried out (cf. [26] for details). For technical reasons we 
apply the Hamiltonian in the first quantized form. The eigen- 
value problem reduces to solving the stationary Schrodinger 
equation H^> = E 1 ^, with H = J2 a H a + Hab composed 
of the single species Hamiltonian 

i<j 

and the inter-species coupling part 

N A N b 

Hab = X! 9ABS <T (xx i - x Bj )- 

<=1 3=1 

Here the effective ID contact interaction potential is mol- 
lified with a Gaussian 5 a (x) = e~ x l 2rT /VZna (of width 
(7 = 0.05) for numerical reasons. In the further examination 
we focus on the case of a harmonic confinement, U(x) = \x 2 
and on attractive inter-species forces qab € (— oo, 0], (The 
case of repulsive inter-species couplings has already been 
investigated in Col . Note that, in the case of U — and 
(la = 9ab, this system is integrable via Bethe's ansatz as in 
1.) 

Computational method Our approach relies on the nu- 
mericallyexact multi-configuration time-dependent Hartree 
method [38,113 EH > a quantum-dynamics approach which has 
bee n ap plied successfully to systems of few identical bosons 
|[ll[Il[l9|,EIlElEl as well as to Bose-Bose mixtures JH]. 
Its principal idea is to solve the time-dependent Schrodinger 
equation i%r(t) = H^(t) as an initial-value problem by ex- 
panding the solution in terms of direct (or Hartree) products 

AN). 
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The unknown single-particle functions <p\ ' (j = 1, 
are in turn represented in a fixed basis of, in our case, 
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harmonic-oscillator orbitals. The specific feature of the 
system at hand is the indistinguishability within each 
species. Therefore the single-particle functions are identi- 
cal within each subset Ka = {1, ...,Na} and Kb = 
{N A + 1,...,N} (, i.e. iff = yf\ Vac € K a ). The per- 
mutation symmetry within each subset is ensured by the cor- 
rect symmetrization of expansion coefficients Aj. In analogy 
to the wave function, also non-separable terms of the Hamil- 
tonian such as the two-body interaction are expanded in terms 
of direct products iRill . 

Note that, in the above expansion, not only the coefficients 
Aj but also the single-particle functions ipj are time depen- 
dent. Using the Dirac-Frenkel variational principle, one can 
derive equations of motion for both Aj, ipj l38ll . Integrating 
this differential-equation system allows us to obtain the time 
evolution of the system via (fj). This has the advantage that the 
basis set {&j(t)} is variationally optimal at each time t; thus 
it can be kept relatively small. Still, its exponential growth 
with the number of particles limits our approach to only a 
few atoms N < 10, depending on how many single-particle 
functions need to be included to describe inter-particle corre- 
lations. 

Although designed for time-dependent studies, it is also 
possible to apply this approach to stationary states. This is 
done via the so-called relaxation method. The key idea is to 
propagate some wave function *(0) by the non-unitary e~ Hr 
(propagation in imaginary time.) As r — > oo, this exponen- 
tially damps out any contribution but that originating from 
the true ground state like e -( E ™-E )T ^ j n p rac tice, one relies 
on a more sophisticated scheme termed improved relaxation, 
which is much more robust especially for excitations. Here 
(\f r |ii"|\f f ) is minimized with respect to both the coefficients 
A j and the orbitals ipj. The effective eigenvalue problems 
thus obtained are then solved iteratively by first solving for 
Aj with fixed orbitals and then 'optimizing' ip j by propagat- 
ing them in imaginary time over a short period. That cycle 
will then be repeated. 



III. MIXTURE OF TWO REPULSIVE COMPONENTS 

In this section we investigate two repulsive components 
(g a > 0) with increasing inter-species attraction gAB € 
(— oo, 0]. We start with components of equal intra-species set- 
tings, such as equal intra-species interaction strengths, gA = 
gs, and particle numbers Na = Nb and discuss subsequently 
the changes in the system's behavior when relaxing these con- 
ditions. 



A. Mixture of two fermionized components 

The starting point is the highly repulsive limit of the compo- 
nents, i.e. two quasi-fermionized states within the two species 
with the inter-species interactions gA = gB = 25.0. For small 
inter-species attraction gAB = —0.001 the system is well de- 
scribed by the uncorrected product of two Tonks-Girardeau 
(TG) states * = ^ A ® where $a = *s « |*o !■ This 



Figure 1: One-body density p{x) = p' a ' (x) (a £ {A, B}) for a 
quasi-fermionized mixture (g„ = 25.0) with the particle numbers 
N a — 2, plotted for different inter-species interactions gAB (see 
legend). 



means, the state of each species a in the high-interaction limit 
(.9a - > +oo) can be mapped to a non-interacting state of 
identical fermions [9], also commonly termed fennionization. 
By extension, a mixture of two fermionized Bose gases has 
similarities with a two component Fermi gas. 

The characteristic fermionic pattern of the TG-state is dis- 
played in the one-body density (which is the same for both 
species a for symmetry reasons) p (x) = p^ (x) (with 
a e {A, B}), which measures the probability distribution of 
finding one a particle at the position x pictured in Fig. Q] 
One recognizes the density concentrating more and more 
in the center of the trap with increasing inter-species attrac- 
tion gAB- Concurrently, in the intermediate interaction regime 
Iff/is | = 5 the initial two density peaks are at first getting in- 
creasingly pronounced, whereas the density in the center of 
the trap grows slowly. That intensification of the fermionic 
characteristic is due to molecule formation, as discussed in 
the following. By contrast, in the very high interaction regime 
(IffAsI = 20 ~ g a ) the two peaks merge into one single peak 
in the center of the trap. 

A more detailed insight in the systems behavior is given by 
the two-body correlation functions. In the case of a binary 
mixture, these are defined as 

Pact ( x 1t x 2) N a (N a — l) 

PAB (X A ,X B ) = NjNz^ A ( X a)^b M^B^b)^ ( x a)) 

The two-body correlation functions p aa (xi, X2) and 
Pab (xa,xb) in Fig. [2] depict the conditional probability 
of measuring the one a-particle at the position x\ and the 
other at the position x<z, and likewise for the different species, 
Pab (xa,xb)- The left column in Fig. |2] indicates that 
the two species keep their "fermionic" character, i.e., the 
probability of finding particles of the same kind at the same 
position {xi = x 2 } stays very low up to high inter-species 
attractions (|<?ab| ^ 20). But according to pab (xa,xb) 
(right column), the two species concentrate more and more 
at the same position, which means on the diagonal in two 
separate peaks aside the center of the harmonic trap (see Fig. 
El Pab (%a, %b) f° r 9ab = —5.0). This can be understood 
as formation of a molecular Tonks-Girardeau (MTG) state: 




/ 




Figure 2: Two-body correlation functions p aa (xi,X2) (left col- 
umn) and pab (xa,xb) (right column) of two quasi-fermionized 
components qa = Qb = 25.0 for inter-species couplings gAB = 
-0.001, -5.0, -20.0, -30.0 (from top to bottom). 



As we shall argue below, two distinguishable particles form 
a bound state, that is a molecule (in the following denoted 
as AB-molecule). These indistinguishable AB-molecules in 
turn form a (molecular) Tonks-Girardeau state. Whereas for 
'two-component Fermi gases this MTG state remains stable 
even in the strongly attractive inter-species attraction regime 
i34ll . this is not the case in a pure bosonic mixture (see also 
Kg. 13. 



Pairing description 

For a better understanding of this behavior we examine the 
Hamiltonian for the exemplary case N a — 2. To this end, we 
transform X = (xa 1 , x A2 ; x Bl , x B2 ) T to the relative coordi- 
nates Y — (Rc m , Ri , J"i j Ti ) T specified by 
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The coordinates i?cM (n,, r 2 ) coincide - up to a factor - with 
the standard center-of-mass (inter-species relative) coordi- 
nates. The coordinate Ri = \ \{xa x +sbJ — (xa 2 + xb 2 )] 
specifies the distance between the centers of mass of two 
(A, B)-clusters. The orthogonal transformation leads to the 
Hamiltonian H (Y) = hcM (Rcm) + H re i, with 
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(4) 

If we assume the formation of, say, AiBi-bound states (i E 
{1, 2}) (up to permutation symmetry), for high enough inter- 
species attraction qab the extension of an A. L Bi- molecule is 
much smaller than the distances between two such molecules 
(| r^l <C |i?i|). One can check this in the two-body correlation 
functions (Fig. [2J. In this limit, we can approximate © by 
the decoupled Hamiltonian 



H, 



re! 
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■gS(Rx)] , 



(5) 



the last part describing the relative motion of the two AB- 
molecules with the effective interaction g = gA+ g_B + 2gAB- 
The analytic solution of the ground state is known [45] and the 
relative part (excluding the trivial CM factor) can be written 
as lOT 




ibrel {X) OC S+ ■( I n fi-^K" 

x U (-e (g) , \ \{x Ax + x Bx ) - {x A2 + x B2 )]) } , 

(6) 

where e(g) = v[g) + 1/2 is determined by the tran- 
scendental equation v (g) £ f^ 1 (0) : fg{ v ) '■= 
2 3/2 [r(i^) /T (-§)]+# and C7 (a, 6) denote the parabolic 
cylinder functions. The symmetry operator S+ := S 1 ;^ <g> 5+ 
serves to compensate the symmetry breaking introduced in the 
Hamiltonian ©. 

This solution gives a good approximation of the den- 
sity patterns in Fig. [2] for intermediate to strong attrac- 
tions (|<?ab| < 3a/2). Also in the high coupling regime 
{\gAB I > <?a/2) the model provides applicable predictions for 
the system's behavior. Considering the molecule-molecule 
interaction term g5{R\) = (gA + g B + ^9ab) 8 (Rx), 
with large enough inter-species attraction gAB the effective 



molecule-molecule interaction g vanishes and even becomes 
negative, i.e. attractive. That implies, more precisely, that for 
9ab ~ 5q/2 a state forms where the AB— molecules are con- 
densed similar to a Bose-Einstein condensate (BEC). For fur- 
ther increase of the interaction, the gas of ^IB-molecules col- 
lapses and forms a bound state. Even though the introduced 
approximation model gives reasonable predictions in that 
limit of very high inter-species attractions (|<7ab| > Sa/2), it 
should be handled with care, as the inter- and intra-species 
length scales become comparable and therefore the scale sep- 
aration breaks down (see also Fig. [2j. 

These considerations can be supported by means of the 
quantities of the one-body density matrix p± (x, x') = 

a {x) (x')) and the pair density matrix 



p(x,x') := 



N A N B 

(A\ B (x)A AB (x')) , 



with the "pair" operator Aab (x) = if? a (x) ^b (x), that 

annihilates an AB-paii "particle" at the position x. As 

Aab (x) \*f?) is a "hole"-state, i.e. a state where an AB-pair 

has been removed at the position x, the pair density matrix 

embodies the overlap of two such "hole"-states. The pair 

density matrix reflects the correlation inherent in the state W 

between the positions x and x' on the level of AB-dimeis, 

as opposed to correlations of single particles a described by 

( Q ) 1 f\ 
Pi (x,x ) ■ 

As shown in Fig. [3] the off-diagonal range of the pair den- 
sity matrix p(x 1 x') persists and even slightly increases for 
the inter-species attractions up to |<?ab| ~ 5, where its ap- 
pearance agrees well with the corresponding one-body den- 
sity matrix /9j [x, x') of identical, fermionized bosons with 
mass M = 2 (in units of to). This proves the existence of 
a paired state (MTG), as discussed above. By contrast, the 
single-particle density matrix p^ (x, x') (see Fig. [3]left col- 
umn) shows two peaks on the diagonal, while the off-diagonal 
density steadily diminishes with increasing |<7ab|- in this light 
a single a-atom will be in an incoherent superposition of left- 
(right-) localized states, without any phase correlations. This 
has to be contrasted with the phase correlations present for the 
pair density matrix (see Fig. [3]right column for gAB = —5.0). 
Interestingly this may be compared to a demixed state in the 
presence of repulsive inter-species interactions 12611 . 

When further increasing the inter-species interaction to 
\gAs\ ~ 20 the size of the system decreases; however, a 
seemingly perfect off-diagonal long range order [46] in the 
pair density matrix is attained, which can be interpreted as 
(few-body analog of) a condensed state on the level of ^IB- 
molecules, p(x,x') — ip* AB (x) ■ (Pab(x') (see Fig. [3]right 
column for gAB = —20.0). By contrast, on the single-particle 
level, displayed in the one-body density matrix p^ (x,x'), 
no condensed state exists, but the two correlation peaks merge 
into one centered peak concentrated on the diagonal {x = x'}. 

For inter-species interaction strength larger than the or- 
der of magnitude of the intra-species interactions \gAB\ ^S> 
\g a \ /2 the system becomes highly bound beyond the AB- 
molecule level, as is reflected in the decrease of the off- 
diagonal elements in the pair density matrix p(x,x'). This 




Figure 3: One-body density matrix p\ (x, x') (left column) and pair 
density matrix p (x, x') (right column) for inter-species couplings 
9ab = -0.001, -5.0, -20.0, -30.0 (from top to bottom). 



can be seen as a collapse from a molecular gas to a strongly 
interacting cluster of AB-molecules. In contrast to a Bose- 
Fermi 13711 and Fermi-Fermi mixture ll3~4ll the collapse in a 
pure bosonic mixture is qualitatively different: In a Bose- 
Fermi mixture only the bosons form a small region with high 
density, whereas the fermions will be attracted up to a "Pauli- 
allowed" density, and Fermi-Fermi mixtures with s-wave in- 
teractions remain mechanically stable even in the strongly at- 
tractive inter-species regime. 

Although, for definiteness, we restricted our discussion to 
a mixture with particle numbers N a = 2, the mechanism 
described by Eq. (0 extends to the case of more molecules, 
N a > 2. For much larger systems N ^S> 1, of course, few- 
body effects as the density oscillations seen in Fig. Q] will be 
smeared out and the corresponding density profiles broadened 
due to repulsion, as in the single-component Bose gas j47h . 
Likewise, the collapse witnessed for qab — > — oo will be 
much more pronounced, as the center-of-mass width Ai?cM 
shrinks with increasing atom number. 



Figure 4: Pair density matrix p (x, x') for g a = 25.0 (left column) 
and the super- Tonks-Girardeau state g a = —15.0 (right column) 
(a 6 {A, B}). The inter-species interaction parameters are qab = 
—0.01, —5.0 (from top to bottom). 



Fermionized attractive components 

It is known for identical bosons that fermionization can also 
be obtained in the attractive interaction regime jl7lll9ll . where 
it is called the super-Tonks-Girardeau state (STG). We show 
that the above pairing mechanism can also be observed in a 
mixture with two attractively interacting, fermionized com- 
ponents (g a < 0, a S {A,B}). In this case it is no longer 
the ground state but an excited state of the system. We per- 
formed the numerical investigation for the exemplary case of 
a mixture with N a = 2 being situated in the energetically 
lowest STG-state (g a = —15.0). Direct comparison with a 
system of repulsive, fermionized components shows the cor- 
responding process analogous to the formation of the molecu- 
lar TG-gas, but with smaller off-diagonal correlations (Fig. [4]). 
Clearly the density profile of the STG-state is more local- 
ized in fragmented regions than that of the TG-state (Fig. |4|. 
The reason is the finite intra-species interaction-strength (here 
g a = —15.0), where the state is not completely fermionized. 
Since this quasi-STG-state still has a non-vanishing, positive 
ID-scattering length a a — — j- > 0, it is closer to a gas of 
spatially extended, hard-core particles (so-called hard rods) 
than to a completely fermionized, point-like TG-gas 11711 . and 
localization effects are more pronounced, which can be ob- 
served in the more profiled density. 



B. Weakly interacting components 

Following the pathway to weak intra-species repulsion, the 
mechanism of pair formation is getting constantly weaker till 
it vanishes in the weak-interaction regime (g a ~ 1). In this 
weakly interacting regime we turn to the limit case of two 
hardly interacting, BEC-like components (g a w 0). Com- 
pared to the case of two strongly repulsive components, there 
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Figure 5: One-body density p (x) of = 2 point-like molecules of 
mass M = 2 with interaction strength g = (dashed), two non- 
interacting components with gAB = —10.0 (solid thin), and two 
molecular components with gAB = —10.0 (solid thick) . 



is no formation of a condensed state of AS-molecules, but 
the system collapses with the increase of the inter-species at- 
traction. In other words, between the AS-molecules there is 
always an effective attractive interaction and thus for strong 
interaction a bright-soliton-like state evolves. Figure [5] dis- 
plays a comparison of the one-body densities of (z) two iden- 
tical point-molecules each of mass M = 2, which mirrors the 
case of very tightly bound, point-like AiJ-molecules with no 
molecular interaction g = 0, (ii) an TV -body bound state of the 
form 

-^sif E K-*B 3 .f) „ ol| . 

$o(-R){e V.is»' 3 J J] e = K"^! } 

a£{A,B} 

and (Hi) the case g a = 0, which is a solitonic state in between 
the two extremes. 

We note that the coherence between the "AS-molecules" 
(as evidenced in the pair density matrix) is slightly stronger 
compared to the one-body level (x,x'), as there is just 
explicit interaction between the species (gAB)- However, as 
there is no longer a scale separation, one cannot consider this 
system simply as a gas of point-like molecules. 



C. Imbalanced components 

After having studied the mechanism for equal component 
settings, we now want to highlight the effects of relaxing the 
equality of the particle numbers and the intra-species interac- 
tion strengths. 
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Figure 6: One-body density matrix p\ ' (x,x') of two quasi- 
fermionized components (g a — 25.0) for the inter-species interac- 
tion strength gAB = -0.01, -10.0, —20.0 (from left to right) of a 
mixture with particle numbers Na = 3 and Nb = 2 (upper row), 
Na = Nb = 3 (lower row). 

Na = 3 and Nb = 2. On the way from weak to very strong 
inter-species attractions gAB, an analogous pathway occurs 
as for fermionized binary mixtures with equal particle num- 
bers, as can be checked on the basis of the pair density matrix. 
That is, a MTG state forms in the intermediate inter-species 
interaction regime, followed by condensation and collapse for 
even higher inter-species attractions. The effect of the dif- 
ference in the particle numbers (or particle densities) can be 
seen as a formation of two phases: One consists of tightly 
bound AS-molecules, as in the case of equal species numbers, 
and the other consists of = \Na — Nb\ (here: Nd = 1) 
"loosely bound" spare particles, i.e. particles that are hardly 
affected by the inter-species interaction (of course, taking into 
account the proper particle exchange symmetries.) This pic- 
ture of loosely bound particles provides a good understanding 
of the two-body density patterns in the intermediate to strong 
inter-species interaction regime (|<7ab| ^ 10). Furthermore 
this formation of, in this case, two j4i?-molecules and one 
loosely bound particle manifests in the one-body density ma- 
trix p^ 1 (x, x') (Fig. [6]) in the formation of two density peaks 
on the diagonal {x' = x} and the larger off-diagonal density 
compared to the balanced counterpart (Na = Nb = 3), re- 
spectively. This two-phase picture breaks down as the system 
starts to collapse for larger attraction (see Fig. [6] gAB = —20). 

Unequal inter-species repulsions 



Unequal particle numbers 

We first consider the case of unequal particle numbers 
Na 7^ Nb, but still the same intra-species interaction 
strengths (gA —9b)- We exemplify this on the case of two 
quasi-fermionized species (g a = 25.0) with particle numbers 



Now we consider unequal intra-species repulsions gA ^ 
gs- For the sake of clarity, we keep the particle numbers 
equal Na = Nb- We start with the case of all intra-species 
interactions corresponding to the fermionization regime, here 
gA = 10.0 and gs = 25.0. In the chosen example the species 
A is weakly fermionized, but still the system evolves simi- 
lar to the case of two strongly fermionized species discussed 
above, that is we observe the formation of a pronounced 
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(c) (d) 

Figure 7: Pair density matrix p(x,x') for the particle numbers Na = 
N B = 2 plotted along the diagonal {a;' = a;} for x £ [-2.2, 0.0) 
and along the off-diagonal {x' — —x} for x 6 [0.0,2.2], of one 
fermionized component gs = 25.0 and intra-species interactions of 
the other component: (a) qa = 0.01 (b) gA = 5.0 (c) gA = 10.0 (d) 
gA = 25.0. The highest inter-species interaction (solid line) shows 
the best possible symmetry between on- and off-diagonal. (The den- 
sities have been rescaled to the same maximal value at the position 
x = 0.) 
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(|<?Ab| > 5), the system collapses in a way characteristic 
for Bose-Fermi mixtures 1 37 ] . That is, the less repulsive 
component A forms a high density region in the center of 
the trap, whereas the strongly fermionized component keeps 
its fermionic character up to higher inter-species interactions 
{\gAB | ~ 8.0). That decrease of the density overlap is visual- 
ized in the light increase of the distance Dist of the extrema of 
the one-body density for 5 < \qab\ ^ 8 (Fig. [§). This char- 
acteristic does not hold for very high inter-species interactions 
(IsabI ^ 10), where the system collapses similarly to that of 
two equally fermionized Bose components, discussed above. 
The comparison of the case at hand (see Fig. |7](c)) with that 
of two fermionized components (see Fig. 0(d)) shows that the 
length scale of the off-diagonal profile is always smaller than 
that of the diagonal profile. Thus no gAB exists for which a 
perfect off-diagonal long range order is achieved in the pair 
density matrix p (x, x'). In this sense, the system starts to col- 
lapse, without forming a condensed state on the AB-molecule 
level. Noteworthy, if both components are equally weakly 
fermionized (gA = gs = 10.0), no partial collapse occurs as 
in the above case, and a condensed state can be observed. 

Going towards components with intermediate and weak 
intra-species repulsions (like gA < 5 shown in Fig. [7] (a- 
b)) the pair-formation is not visible anymore, and the system 
immediately starts to collapse without forming a condensed 
state as in the case before. In the extreme case of one "con- 
densed" component (17^4 — > + ), the approximation (O above 
is not valid as the length scales cannot be separated anymore. 
For unequal particle numbers, the system can again be thought 
of as a two-phase system; hence if the particle number of the 
fermionized state (say B) exceeds the number of condensed 

( B) 

particles, the coherence in the one-body density p\ (x, x') 
has larger off-diagonal elements due to spare (unbound) B- 
particles. 



I9abI 

Figure 8: Distance Dist between the one-body density maxima 
Max^"' (x, x')] of the components a 6 {A, B} for different inter- 
species interaction strength gAB and the intra-species interactions 
g A = 10.0 and g B = 25.0. 

Tonks-Girardeau pattern in the pair density matrix p (x, x'), 
indicating the MTG state. This is in line with Figure|7l which 
pictures the p (x, a;') -profile along the diagonal {x' = x} and 
along the off-diagonal {x' = — x}. 

During the formation of the MTG state, one can observe 
an assimilation of, for instance, the one-body density matri- 
ces p^ (x, x') and p[ B ^ (x, x'). The best overlap is achieved 
about the value |<?ab| ~ 5.0. To characterize this increas- 
ing similarity, let us define the distance between the peak po- 
sitions x max,a of each component's (diagonal) density pro- 
file, i.e., where p± (x( max ' a \ j( max i a )) is maximal). For 
for a given inter-species coupling gAB, this is denoted by 
Dist := || x (""^)| - ^(max,^ | ( see Fig. [H>. 

With a further increase of the inter-species interaction 



IV. MIXTURE WITH ATTRACTIVE COMPONENTS 

In this section we complete our investigation by exploring 
mixtures with one or more attractively interacting components 

(g a < 0). 



A. Repulsive and attractive components 

We start in the spirit of the above section with one compo- 
nent in the fermionized interaction limit, i.e. gs = 25.0 and 
the other in a bound state gA = —10.0. The bound species are 
strongly localized in the center of the trap and the feedback on 
that component is negligible for any inter-species attraction. 
To explore this situation it is natural to consider a simplified 
Hamiltonian, where the effect of the localized species A is 
replaced by an additional external potential for the B atoms 
8U B {x) = g A BN A p (A) (x) JH. Here we apply the analyti- 
cally well studied split-trap p^ (x) w <5 (x) 11481 l49ll . 
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x 1 x 1 x 1 

Figure 9: Two-body correlation function pbb (asi, X2) of a mixture 
with one molecular species qa = —10.0 (Na = 4) and one species 
of Nb = 2 repulsive bosons, gs = 25.0, for the inter-species inter- 
action parameter qab = —0.01, —1.0, —5.0 (from left to right). 
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Figure 10: One-body density p' s ' (x) of Nb = 2 non-interacting 
B-particles (gs = 0.0), with different A-particle numbers Na in a 
bound state (gA = —10.0) and the inter-species interaction strength 

JAB = -10.0. 



h b = E {\p 2 b + \ x b + N A g AB S (x Bi )) 

(7) 

+9b E $ (XBi - x Bj ) ■ 

Furthermore for gs ^> 1, one can map the fermionized 
component on a non-interacting fermionic system [9]. Con- 
sequently one obtains for the exemplary case of Nb — 
2 S-particles the simple solution 1^ = |^rro«m| = 

-^\®o{xx)$x{x 2 )-®o (x 2 )^i(xi) |, where $, denotes 
the i-th single-particle eigenstate of a split-trap. 

The validity of the approximation becomes rapidly better 
with increasing number of particles (Na) in the bound state, 
as the width of one-body density (x) scales as ^== , and 
hence converges towards a <5-type potential in the limit of large 
particle numbers (Na — > 00). The agreement is astonishingly 
good already with relatively few particles Na > 4. In Fig. 
[9] the exact two-body correlation function pbb (xx,X2) for 
Nb = 2 particles is shown. As it turns out, the picture of a 
non-interacting fermionic system applies very well up to in- 
termediate inter-species interactions \gAs\ < 2 (in the case of 
Na = 4). This predicts that if for intermediate inter-species 
interaction qab ~ — 1 one detects a £> -particle aside the trap 
center (xb x ~ ±1), the other £>-particle is located at the cen- 
ter of the trap (xb 2 ~ 0). (For Nb > 2 an additional density 
contribution emerges on the off-diagonal.) 

Whereas for the model (5-type potential the two-body den- 
sity pbb (xi,X2) would remain in the (increasingly sharp) 
cross-shaped pattern even for qab — ► —00, this is not the 
case for the system at hand (see Fig. [9] for qab = —5.0). 
Due to the nonzero width of the additional potential caused 
by the A-particles, for high enough inter-species attraction all 
of the B particles can be accommodated in the "A-potential" 
as a whole (unlike for a (5-type potential). That is illustrated in 
Fig. |9]for qab — —5.0, where again the familiar TG-density 
pattern can be observed, but with the spatial extension of the 
y4-particle density (compare with Fig. ITOb . This behavior in 
the high interaction regime can also be observed for higher 
B-particle numbers. 



B. Two attractive components 

Extending the results of the last section, we start with the 
case of one weakly interacting, i.e., "condensed" component 
(c)b ~ 0), and the other component again in a bound state 
(gA = —10.0). We can again apply the previous split-trap ap- 
proximation on a system with Nb = 2, i.e. the condensed 
S-particle feel an effective short-range potential at the center 
of the harmonic trap. Again with increasing particle-number 
Na in the molecular state, the approximation is getting bet- 
ter. However, for the length scales of the two components to 
differ as distinctly as in the case of a fermionized component, 
the agreement with the split-trap model above requires more 
S-particles in the bound state. If we assume as a model a con- 
densed state in component B (gB = 0) and a tightly bound, 6- 
type state in component A, the model Hamiltonian ^} reduces 
to 

N B 

Hb = E (b 2 Bi + h x %+NAg A BHxB t )) , (8) 

i=l 

with the solution JH] * (X B ) oc 
ex P (-^4^)nf=t U (i - f , lx%), which 

evolves with increasing inter-species attraction towards a 
state analogous to the bound state of the ^-potential. 

For higher attractive interactions also in component B, the 
model is not applicable anymore. In the limit of highly at- 
tractive components, with increasing inter-species interaction 
strength the system forms an entire bound state. For the spe- 
cial case of equal interaction strength g a = gAB = g, one can 
map the system to a bound state of Na+Nb identical particles 

*o (X = (X A , X b )) cx $ (R) exp ( - M £ \ Xl - Xj \ j . 



V. SUMMARY 

In conclusion, we have investigated the ground state 
of a two-component Bose gas in a one-dimensional har- 
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monic trap throughout the crossover from weak to strong 
inter-component attraction. We have highlighted different 
pathways depending on the choices of the different intra- 
component couplings, and indicated how they can be un- 
derstood in terms of simplified models. For two quasi- 
fermionized components (i.e. Tonks-Girardeau states), the 
system forms a molecular Tonks-Girardeau gas in the inter- 
mediate inter-component interaction regime, which consists 
of bound pairs containing one particle of each component. In 
the strongly attractive regime, we demonstrated the condensa- 
tion of the bound pairs, followed by the collapse of the system 
beyond a critical attraction. We showed the analogous mecha- 
nism for attractively fermionized components, that is compo- 
nents in the super- Tonks-Girardeau regime. Relaxation of the 
condition of two equally fermionized components leads to a 
modified pathway: In the case of just one fermionized com- 
ponent, the formation of a molecular Tonks-Girardeau gas can 
still be observed for high enough repulsion within the second 
component, but the collapse occurs in analogy to Bose-Fermi 
mixtures without condensation of the bound pairs, in contrast 
to the case of comparable repulsions. In the regime of inter- 



mediate (inter-species) attraction, unequal (number-) densities 
in the components have been found. These can be understood 
as two phases, one consisting of molecular pairs of each com- 
ponent, and the other phase consisting of loosely bound par- 
ticles. For mixtures with one strongly attractive component, 
we showed that this component can be interpreted as an addi- 
tional external <5-function potential for the other component, 
in case both length scales can be well separated. The investi- 
gation of these intriguing pairing scenarios paves the way to- 
ward studying their quantum dynamics, such as the tunneling 
of molecular pairs in multi-well traps. 
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